Quantum dimer phases in a frustrated spin ladder: Effective field theory approach 

and exact diagonalization 
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The phase diagram of a frustrated 5=1/2 antiferromagnetic spin ladder with additional next- 
nearest neighbor exchanges, both diagonal and inchain, is studied by a weak-coupling effective 
field theory approach combined with exact diagonalization for finite systems. In addition to two 
known phases with rung-singlet and Haldane-type ground states, we observe two new phases with 
dimerization along the chains. Furthermore, the transitions between the different phases are studied 
and shown to be either first order or to belong to the universality class of the two-dimensional Ising 
model. The nature of elementary excitations is discussed briefly. 
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PACS numbers: 75.10.Pq, 75.40.Cx, 75.30.Kz, 75.40.Mg 



INTRODUCTION 



During the past few decades, strongly correlated elec- 
tron systems have received a lot of attention, e.g. due to 
unconventional superconductors. In particular, during 
the early days of high-T c superconductivity, Anderson 
had already proposed a mechanism based on a so-called 
'resonating valence-bond' (RVB) picturei Nevertheless, 
the search for a spin liquid in two dimensions remains 
a long standing problem. In this context, much atten- 
tion was devoted to the study of geometrically frustrated 
magnetic systems (for recent reviews see e.g. Refs.00)- 

Among the different proposed models to capture An- 
derson's RVB scenario is the one of Nersesyan and 
TsvelikjA based on a spatially anisotropic Ji — J2 square 
lattice with stronger exchanges along one particular 
'chain' direction to allow for generalizing results obtained 
previously in one dimension. 5 However, it turned out that 
even the one-dimensional picture was not completely un- 
der control^ in weak coupling it is impossible to fine-tune 
rung and diagonal exchanges of two-leg spin ladders to 
eliminate all relevant interchain couplings at all orders. 
In particular it was argued in Ref. |fj that at weak cou- 
pling there is no direct first-order phase transition from 
two phases with unique ground states, namely a Haldane 
phase and a rung-singlet phase, rather an intermediate, 
spontaneously dimerized phase was predicted. However, 
previous numerical calculations have not detected this 
intermediate phased The numerically determined phase 
diagram suggested a direct phase transition from the Hal- 
dane to the rung-singlet phase which at weak interchain 
coupling seemed to be of second order and become first 
order at stronger couplings. While the latter topology of 
the phase diagram agrees with the earlier bosonization 
prediction) 5 " the order of the phase transitions disagrees 
at weak coupling. We will argue that the intermediate 
phase in the spin ladder can be revealed by adding fur- 
ther exchanges. 



In this work we study an S = 1/2 antiferromagnetic 
two-leg spin ladder with additional next-nearest neighbor 
exchanges and discuss the interplay between two distinct 
scales, one set by inchain frustration and the other one by 
interchain interactions. In the spirit of Ref.|8, this model 
can also be regarded as a strip of the extensively studied 
J\ — Ji — J3 square latticeAlft The latter model illustrates 
that frustrated interactions in quantum spin systems give 
rise to 'unpredictable' phases which are very difficult to 
analyze theoretically because of the presence of compet- 
ing interactions. The study of the effect of additional 
next-nearest neighbor exchanges in spin ladder systems 
therefore also is an interesting theoretical problem in its 
own right. The bosonization approach 1 ^ is an unbiased 
and powerful method for studying frustrated systems in 
particular in one dimension. Using this method, we will 
show that a rich phase diagram and intriguing nature of 
elementary excitations result from the competing inter- 
actions in the frustrated spin ladder. Furthermore, nu- 
merical analysis of finite systems confirms the presence 
of the phases obtained within the weak-coupling effective 
field theory, also at larger couplings. 



II. MODEL 

We consider the S — 1/2 antiferromagnetic Heisen- 
berg spin ladder with additional next-nearest neighbor 
exchanges which are represented by ladder diagonal and 
inchain next-nearest neighbor interactions. The geome- 
try of our model is depicted in Fig. ^ The lattice Hamil- 
tonian reads: 



= Hleg + Hleg + Hint , 

where the Hamiltonian for leg I = 1 , 2 is 
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FIG. 1: Structure of the spin ladder with next-nearest neigh- 
bor interactions. 



and the interleg coupling contains both the rung and the 
diagonal exchanges: 



JY 



Hint — [J±. SlJ • S2j" 

+Jx (Si j- ■ S2J+1 + Si J+ i • S2J; 



(3) 



Above Si j represent spin 5=1/2 operators on the j- 
th rung, l-th leg, and periodic boundary conditions are 
assumed S/^at+i = Sj,i. In order to avoid additional frus- 
tration by the boundary conditions, we will consider only 
even chain lengths N. 

The above Hamiltonian Q with J2 = has already 
been investigated in detaili£ii2iiii In particular the case 
J x = J, J2 = serves as a good illustration of two 
possible phases^ since in this case the total spins on all 
rungs are good quantum numbers. Then the ground state 
consists of local singlets S run g = for sufficiently large 
J_l with a clear gap to all excitations. On the other hand, 
for small J± , total spin SV ung = 1 is formed on each rung 
and the low-energy physics is governed by an effective 
spin-one chain. In the latter case, the ground state is the 
famous 'Haldane' state which is also gapped (see e.g. Rcf. 
fill for a recent review of spin-chain models). These two 
phases extend to J x ^ J, J2 7^ and we refer to them 
as 'rung-singlet' and 'Haldane' phase, respectively. 



III. WEAK-COUPLING APPROACH, 
BOSONIZATION 



In this section we perform the weak-coupling analysis 
of our system (I) . We will follow the usual way and start 
from a continuous field theory description of the indi- 
vidual spin chain and treat the interchain interactions 
perturbatively^i We note that one can construct two 
equivalent weak-coupling formulations. One can start 
from decoupled chain limits where chains run either along 
the ladder legs, or along the ladder diagonals. These two 
cases are connected to each other by a 'duality trans- 
formation' exchanging the spins on every second rung^ 
For definiteness we start from two decoupled frustrated 
chains running along the ladder legs: J x , J± <C J. Spin 
operators on each chain are decomposed in their smooth 



and staggered parts: 

S,(x = jao) = J t (x) + 3i(x) + (-lfn,(x) (4) 



aa 

with lattice constant aQ. J; and J; are left and right 
SU(2) currents of the l-th chain and have conformal 
weights (1,0) and (0,1) respectively. Staggered spin den- 
sities are represented by more relevant operators, indicat- 
ing an inherent instability towards doubling of the unit 
cell in antiferromagnetic chains: 



n\{x)~Tr{o l 9l {x)) 



(5) 



Here gi(x) stands for the basic 2x2 matrix field 
of the Wess-Zumino model and has conformal weights 
(1/4,1/4). a 1 are the Pauli matrices. The scalar quan- 
tity which is represented in the continuum limit by an 
operator of the same dimension as the staggered part of 
the magnetization is the dimerization operator: 

Si(x)Si(x+ao) = (—l) x ei(x) + \css relevant smooth part. 

(6) 

In terms of the basic Wess-Zumino matrix field the fol- 
lowing representation holds for the dimerization opera- 
tor: 



ci(x) - Tig^x) 



(7) 



The effective quantum field theory in the continuum 
limit of a single chain is that of a critical SU(2)i Wess- 
Zumino model perturbed by a marginal current-current 
interaction>i£iiii 

fl- I = ^-(:J,J I : + :J,J I :)+7J,J I (8) 



with the following notations: 



and 



7 = J-2 - J: 



2,c : 



where J2, c — 0.24iiLi£ If 7 < the perturbation is 
marginally irrelevant, but when 7 > the interaction 
flows towards strong coupling and the system dimerizes 
spontaneously with a dynamically generated gap in the 
excitation spectrum. 19 In Abelian bosonization represen- 
tation we can rewrite © as: 



Hi 



2vr 2 



[d x <j>i{x)Y 



7 

+ -d x (j)L R d x (f>i L 



(9) 



with compactified dual bosonic fields 9 and (j>. The fol- 
lowing representations hold for the staggered spin density 
and dimerization operatorssii 



ri; (x) ~ [cos v (x) , sin v27r#i(x), sin V27r0;(x) 

£;(x) ~ COS \p2/K(\)i (x) . (10) 
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For further analyses, treating interchain coupling pertur- 
batively, it will be convenient to pass to symmetric and 
antisymmetric combinations of the bosonic fields 

4> ± {x)= Mx)± r- Mx) (11) 
V2 

and define the continuum limit expressions of order pa- 
rameters for columnar and staggered dimerizations of the 
two chains: 

e+ = ei + 62 ~ cos \A''/ > + cos V^4>- > 

e_ = t\ — t2 ~ sin \/Tt(p+ sin . (12) 

The Hamiltonian in the symmetric and antisymmetric 
sectors will contain a marginally relevant intersector cou- 
pling which promotes dimerization within the chains: 

-d x (j) + , R d x (j) +:L + -d x (f>-,Bdx4'-,L 

— — ^ cos\/47r0+ cos\/47r</)_ . (13) 

The intersector interactions ljl3|l show that the vac- 
uum configurations of the symmetric and antisymmetric 
sectors are degenerate: 

a : ^cosa/47t0+^) = ^cosa/47t<?!>-^ = 1 

b: ^cos V47r0 + ^> = ^cos \/47r0_\ = -1 . (14) 

This includes the two-leg ladder dimerized in different 
patterns, namely a: with long-range ordered columnar b: 
staggered dimerizations, respectively. 

The expressions of spin operators in terms of the 
bosonic fields are valid as long as the dimerization is 
very weak. We can then add interchain interactions to 
the frustrated decoupled chains. We classify interchain 
couplings according to their scaling dimension: 

Hint ~ (J± + 2 J x ) (Ji(x) + Ji(x)) (J 2 (s) + J 2 (x)) 
+ (J x -2J x )n 1 {x)n 2 {x). (15) 

The product of the smooth parts of the spin operators is 
translated into marginal operators, whereas the product 
of the staggered parts is represented by relevant operators 
in the effective field theory. 



A. Mean-field separation 

As a first approximation and for 7 = 0, Jj_ — 2 J x 7^ 0, 
one can retain only relevant terms that stem from the 
product of the Neel components of the spin operators. 
In this case the effective field-theoretic Hamiltonian sep- 
arates into two commuting parts in the symmetric and 



antisymmetric basis i2£L2i 

H = H+ +H- 

H + = ^[(d x 6 + f + (d x <P + ) 2 ]-^4-co S V^Mx) 
2 2ir z 

H~ = ^[{d x 6.f + {d x ^{x)f\ (16) 

J±c 2 J±c 2 
+ cos \/47r0- (x) H 5— cos V4tt9- (x) , 



where Ji = J± — 2 J y ■ , c stands for a non-universal nu- 
merical constant and 



± = u (i±JL) 

y 2vr J J 



in weak coupling. 

Now we add 7 > and look at the limit J± , J x <C 7 <C 
J, which allows to assume stability of the massive dimer- 
ized phases of the individual chains against infinetisimal 
interchain perturbation. For Ji = J x = the chains 
dimerize in columnar and staggered patterns which cor- 
responds to pinning of the symmetric and antisymmet- 
ric fields in degenerate vacua with non-zero averages of 
cos y/^n(j)± . One may assume that these averages will 
be non-zero in a finite region of the parameter space af- 
ter inclusion of interchain coupling such that we may 
substitute finite averages for cos ^f^K(f>± in i|13|) . In the 
antisymmetric sector the Hamiltonian density contains 
a relevant cosine of the bosonic field 0_ as well as its 
dual 6L, see (|16f) . By contrast, the dual field cannot ap- 
pear in the symmetric sector for symmetry reasons. The 
classical vacuum configuration of the symmetric field is 
pinned unambiguously according to the sign of interchain 
coupling and thus any small interchain coupling (which 
is relevant) removes immediately the degeneracy of the 
ground state of decoupled dimerized chains (|14fl . The 
picture is as follows: addition of infinitesimally small rel- 
evant interchain coupling confines the massive spinons of 
the individual chains into magnons. At the same time, 
the ground-state degeneracy is lifted and one of the long- 
range ordered dimerization patterns (columnar or stag- 
gered) is selected. Once the symmetric field 4> + is pinned, 
we can perform a mean-field decoupling of the interac- 
tion term stemming from 7. Note that this procedure 
does not depend on the sign of J± — 2 J x 7^ 0, the cru- 
cial assumption is that a relevant interchain coupling is 
present, which provides a confining potential to spinons 
of the individual dimerized chains and selects one of the 
two long-range ordered dimerization patterns. 

At the mean-field level we are left with the following 
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Hamiltonian: 



form in Major ana fermions: 21 



H 



MF = 



H H — cos V4tt0- (x) 



+ ^ ^cos y/4rt<f>-(x)^j cos VAtt4) + (x) 



(17) 



where 7?o stands for the sum of the Gaussian parts of 
the symmetric and antisymmetric sectors. This particu- 
lar form of the mean-field Hamiltonian is convenient for 
passing to Majorana fermions, but let us first discuss ar- 
tifacts of the mean-field separation: the mean-field sepa- 
ration in i|17|) breaks the underlying SU (2) symmetry of 
the ladder model to t/(l)<g)Z 2 . This is verified directly by 
noting that the very form of the mean-field Hamiltonian 
(|17H coincides with the Hamiltonian of a ladder where 
the next-nearest inchain frustrating coupling is substi- 
tuted by the product of dimerization operators of the 
two chains and interchain coupling is of easy-plane type: 



with 



(A + B)e 1 {x)e 2 {x) + {A- B)n\n\ 



A = --t- (cosv / 47T(/) + (a;)) , 



B 



(18) 



— ^cos VAtt^^x)^ 



While the product of dimerization operators respects 
SU(2) symmetry, the appearance of an effective 
anisotropy in interchain coupling is clearly an artifact of 
the mean-field separation and should not be taken phys- 
ically (e.g. splitting of the triplet of Majorana fermions 
into doublet and singlet). With increasing interchain ex- 
change in l|17fl we observe a phase transition represented 
by a self-dual point in the antisymmetric sector: 



J ± c 2 
2tt 2 



7 



—J ( COS V (x) 



J j 



(19) 



Since our approach is a weak coupling one (starting 
from decoupled chains) we are dealing with a fixed-point 
Hamiltonian which is described in terms of two copies of 
the SU(2) 1 Wess-Zumino-Witten (WZW) model that is 
equivalent to an SO (4) WZW model. The 50(4) WZW 
model admits a representation of its generators in terms 
of four real Majorana fermions. 2 - To clarify the symme- 
tries of the stabilized dimerization pattern and to unveil 
the nature of the self-dual point in the antisymmetric sec- 
tor we pass to these Majorana fermions and subsequently 
to Ising variables. 

The mean-field Hamiltonian (|17|) takes the following 



H 



dx 



IV 



— {PRdxPR - PhdxPh) + im s p R p L 



E / dx 



IV i 



(20) 



where 



m t ~ (J± - 2 J x ) +sign(Ji - 2 J x )7, 

m s » -3(J ± -2 J x )+sign(J ± -2J x )7 (21) 

and 7 ~ exp (— 2ir/j) > in weak coupling. The ap- 
pearance of the sign function in the expressions for the 
masses of the Majorana fermions is the most important 
result of the mean-field treatment. Eq. i|21[l shows that 
only the singlet among the Majorana fermions can soften 
by varying interchain interactions. Note that the triplet 
and singlet Majorana fermions can be interpreted as fol- 
lows in the limit of strong rung coupling: the triplet ip a 
corresponds to a rung magnon (a triplet excitation on 
a rung), while the singlet p corresponds to the singlet 
bound state of two magnons. Indeed, we will argue soon 
that in the context of the effective field theory the sin- 
glet Majorana fermion softens at Ising phase transitions. 
On the other hand, we will show in section HV CI that the 
same Ising phase transitions can be traced to the soften- 
ing of a singlet bound state of two magnons in our exact 
diagonalization data. 

It is well-known that a Majorana fermion describes 
the long-distance properties of the two-dimensional Ising 
models with the mass of the fermion proportional to 
the deviation from criticality m ~ (T — T c )/T c , mean- 
ing in particular that positive mass corresponds to the 
disordered phase of the Ising model. The following local 
representations hold in Ising variables for the order pa- 
rameters of the columnar and staggered dimerized phases 



CTiCT2CT3fT4 , 



(22) 



where Ui and cr, are disorder and order variables of 
the i-th replica of the two-dimensional Ising model 
respectively^ 

From the above formulas we read off the following pic- 
ture: if we add to decoupled dimerized chains a relevant 
coupling with positive relevant exchange, J± — 2 J x > 0, 
then four copies of the Ising model will be in the dis- 
ordered state, meaning that the ground state will have 
columnar dimer order. In the case of Jj_ — 2 J x < 0, a 
staggered dimerized ground state will be selected. This 
analysis is supported by the following topological argu- 
ments: the ground state of the ordinary antiferromag- 
netic ladder is parity symmetric, (e±) = 0. A typi- 
cal configuration in the resonating valence bond state of 
the antiferromagnetic ladder is depicted in Fig. [2 with 
rectangular boxes indicating singlets formed between the 
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FIG. 2: Typical configuration in the RVB state of the anti- 
ferromagnetic ladder. 



nearest-neighbor spins and the zigzag line representing 
an effective S = 1 spin formed across the ladder diago- 
nal. 

One can introduce even and odd topol ogic al ('string') 
order parameters, as described in Refs. l2 r j|2(i These 
string order parameters have a simple geometric interpre- 
tation: they count the number of valence bonds crossed 
by an arbitrary vertical line such that the even string or- 
der parameter is non-zero if the number of crossed bonds 
is even and the odd one is non-zero if this number is 
odd. A nice feature of the two-leg ladder systems is 
that short-range valence-bond configurations have defi- 
nite parity string order parameter^ 

In the rung-singlet phase the even string order param- 
eter is non-zero. Since 7 does not introduce coupling be- 
tween the chains both patterns of dimerization (columnar 
and staggered) are energetically equally favorable at this 
level, and from Fig- EI it is clear that the columnar dimer 
configuration will be long-range ordered after adding rel- 
evant interchain coupling with J± — 2J X > 0. Similar ar- 
guments apply to the case J± — 2J X < 0. From the topol- 
ogy of typical configurations of the RVB state one can 
see in this case that singlets formed along the ladder di- 
agonals coexist with staggered dimer configurations iSSiSi 
The string order parameter is odd and upon increasing 
inchain frustration the staggered dimer phase will be se- 
lected for the ground state. Crossing the surface which 
at weak couplings reads J± — 2 J x ~ 0, the string or- 
der parameter jumps from even to odd by a first-order 
phase transition. These findings are summarized by the 
qualitative weak-coupling phase diagram Fig. 

Eq. (|21|l also shows that gaps open linearly as one de- 
viates from the Ising transition lines either in the rung 
singlet or in the Haldane phase directions. 

Let us summarize the mean-field scenario presented 
above and the supporting topological arguments: when 
two dimerized chains are coupled in such a way that rel- 
evant interchain coupling is present (in the model which 
we consider it is always present at weak couplings 6 ) 
it provides a confining potential between the gapped 
spinons of the individual chains and the four-fold degen- 
eracy of the ground state is reduced to two fold. But 
what can we say when interchain coupling is itself frus- 
trated, i.e. when no relevant terms stem from interchain 
exchange and no confining potential is supplied for decon- 
fined gapped spinons (for example, if we add additional 
four-spin interactions)? We answer these issues in the 
next section where we derive a new quantum phase tran- 
sition when both competing interactions (inchain as well 




J x /J 



FIG. 3: Topology of the phase diagram as captured within the 
weak-coupling approach. Phases are as follows: A- rung sin- 
glet, B- Haldane, C- columnar dimer and D- staggered dimer. 
The dashed line represents a first-order phase transition line, 
whereas continuous lines stand for Ising transition lines. De- 
viations of the dashed line at weak coupling from the straight 
line J± — 2 J x and enlargement of the columnar dimer phase 
are due to second order-processes as found in Ref. [g. 



as interchain) are marginal. 



B. Mapping to the Gross-Neveu model 

In this section we fine-tune interchain exchange to 
cancel relevant contributions. As recently discussed by 
Starykh et al.fL in the simple model interchain ex- 
change cannot be fine-tuned to yield only marginal inter- 
actions unless we add four-spin counter terms ^ In real 
systems small four-spin interactions are always present 
e.g. represented by the cyclic ring exchange. So we can 
assume that we start from the line Jj_ — 2 J x =0 and rel- 
evant terms generated in second order are compensated 
with small four-spin counter terms. In that case the low- 
energy degrees of freedom of our fine-tuned model will 
be described by the triplet and singlet massless Majo- 
rana fermions coupled by marginal four-Fermi interac- 
tions alone: 



= Ja d.T[(A + 7 )(^iV^ 



+ ^l^l + ^l^l) (23) 



where A = J± + 2 J x = 2 J±. As long as 7 < A, including 
negative 7 we can repeat the analysis by Allen et al^ 
where neglecting the anisotropy a mapping to the 0(4) 
Gross-Neveu models was used. Based on the one-loop 
RG equations one can check that negative 7 flows to zero 
at low energies and neglecting the anisotropy will not af- 
fect the ground-state properties of the model, neither will 
the excitation spectrum be modified qualitatively. With 
increasing 7 to zero the mapping to the 0(4) Gross-Neveu 
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model becomes only a better approximation and exactly 
when 7 = the 0(4) Gross-Neveu model is the exact low- 
energy effective field theory describing chains coupled by 
frustrated interchain interaction. Further increasing in- 
chain frustration, telescoping 112MI) shows that at the point 
where 

7 = A = 2 Jl (24) 

a new quantum phase transition takes place where the 
zeroth Majorana fermion decouples from the rest of the 
triplet. Low-energy excitations at this quantum critical 
point are governed by a free massless Majorana fermion 
plus an 0(3) Gross-Neveu model with dynamically gener- 
ated mass for the remaining triplet of Majorana fermions. 
In the extreme limit 7 » ) we can neglect the small 
anisotropy between triplet and singlet Majorana fermions 
and recover the 0(4)-invariant Gross-Neveu model, with 
dynamically generated mass. The sign of the mass de- 
fines the dimerization pattern of the model. Elementary 
excitations are deconfined massive spinons which interpo- 
late between the four-fold degenerate vacua of the lattice 
model. This theory describes a first-order phase tran- 
sition line between the columnar dimer and staggered 
dimer phases. 



IV. NUMERICAL RESULTS 

In order to test numerically the predictions for the 
phase diagram and the nature of phase transitions we 
start from the Majumdar- Ghosh point J2 = J/2 in the 
decoupled chains^ Each of the decoupled chains then 
has two degenerate dimerized ground states with a gap 
to the excitations. This behavior is generic in the vicinity 
of the Majumdar- Ghosh point, but the correlation length 
is minimal exactly at J2 = J/2. 31 Such a short correlation 
length is advantageous for numerical calculations since it 
minimizes finite-size effects. 

For two weakly coupled dimerized chains, there will 
be four states which are well separated from higher ex- 
citations. We then need to clarify which ground state is 
formed depending on the coupling between the chains. 
In this context, it is interesting to note that for the 
particular choice Jj_ = 2 J x of interchain coupling at 
J2 = J/2, the two staggered dimer arrangements remain 
exact eigenstates. Indeed, it is straightforward to check 
that interchain coupling terms cancel on the staggered 
dimer states. On the other hand, columnar dimer states 
cease to be exact eigenstates of the full ladder Hamilto- 
nian, thus they can lower their energy by fluctuations. 
Therefore, we expect that two chains with J2 = J/2 
weakly coupled by J± = 2 J x will have a columnar dimer 
pattern in the ground state. 

To determine the ground state phase diagram, we 
have performed exact diagonalization of ladders with up 
to N = 16 rungs (32 spins) using a standard Lanc- 
zos method, exploiting the conservation of total S z and 
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FIG. 4: (Color online) Numerical results for the lowest ex- 
cited states with S = and k = tt at J2 = J/2, J± = 0.2 J as 
a function of J x . Lines are for N — 16 rungs, symbols show 
extrapolations to the thermodynamic limit. The big filled tri- 
angles denote estimates for the locations of the different phase 
transitions. For details compare the text. 



spatial symmetries. We impose periodic boundary con- 
ditions along the chains in order to have translational 
symmetry. Furthermore, there is an exchange symmetry 
of the chains. Important information about the ground 
state is provided by the associated quantum numbers, in 
particular momentum k for translation along the chains, 
and k y for the 'momentum' corresponding to an exchange 
of the chains. On a finite system, the two staggered and 
columnar dimer states both combine to a ground state 
with k = 0, k y = 0. However, the value of k y in the 
first excited state distinguishes between the two cases: a 
combination of staggered dimer states gives rise to k = tt, 
k y = tt, while one obtains k — tt, k y = from the colum- 
nar dimer states. 

Extrapolations to the thermodynamic limit N — > 00 
will be performed using the Vandcn-Broeck-Schwartz 
(VBS) algorithmic^ (with the parameter a = — 1 of 
the algorithm). 



A. Low-lying levels 

To illustrate the generic behavior, we will first discuss 
a cut through the phase diagram at J± = 0.2 J. The lines 
in Fig.0|show the excitation energy relative to the ground 
state for the lowest levels with k = tt, k y = and tt on 
a finite system with N = 16 rungs. On the left side of 
Fig. 01 ah excitations are gapped, consistent with a rung- 
singlet phase. For < J x < 0.18 J, the k — tt, k y = 
level approaches the ground state, as is expected in a 
columnar dimer phase. Note that the point J± = 0.2 J, 
J x = 0.1 J lies inside the columnar dimer phase, as is 
expected since it belongs to the line J± = 2 J x discussed 
above. On the other hand, for J x w 0.2 J, it is the level 
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with k — 7T, k y — 7r which is approaching the ground 
state. This behavior is consistent with a staggered dimcr 
phase around J x ~ 0.2 J. Finally, on the right side of 
Fig. 01 all excitations are again gapped, consistent with 
a Haldane phase. 

The locations of the phase transitions shown by the 
three big filled triangles in Fig. 01 are determined as fol- 
lows. For intermediate J x , the levels shown in Fig. 0] arc 
the only low-lying excitations, implying a direct first- 
order transition from the columnar dimer to the stag- 
gered dimer phase. We take the crossing of these two 
excited levels on a ladder with TV rungs as an estimate 
for the transition point. We then extrapolate the crossing 
points with 6 < TV < 16 rungs using the VBS algorithm. 
This extrapolation leads to the middle triangle in Fig. 01 
Since the dependence on TV is only weak, this extrap- 
olated point is very close to the finite-size estimate for 
TV = 16. 

The transitions to the rung-singlet phase and the Hal- 
dane phase at the left and right boundaries are expected 
to be second-order phase transitions. Such transitions 
can be accurately estimated using the phcnomenological 
renormalization group method*^ The critical couplings 
{J} are determined by the condition that the levels for 
two different system sizes Ni, N 2 satisfy 

JVi AE k= ^ ky {N u {J}) = N 2 AE k= ^ ky (TV 2 , {J}) . (25) 



Rung singlet 



Here 



one should choose k y = and tt for the transitions 



into the columnar and staggered dimer phases, respec- 
tively. Application of (J2SJ) to N± = 14 and N 2 = 16 
yields the two remaining filled triangles in Fig. 0] (we ne- 
glect the dependence on TVi and TV2 since in this case it 
is only weak). 

Finally, the open circles and triangles in Fig.^show an 
extrapolation to the thermodynamic limit of the k y = 
and the k y — tt levels, respectively. This extrapolation 
has been performed with system sizes starting at TV = 4 
rungs using the VBS algorithm. Errors of the extrap- 
olation are difficult to estimate, but should not exceed 
the size of the symbols. In some cases, no convergence 
can be seen in the finite-size data for the higher-lying 
level such that no extrapolation can be performed for 
that level. However, in these regions of J x , the energy of 
the higher level shows a tendency to increase with TV for 
the accessible system sizes. Thus, the higher-lying level 
can safely be assumed to stay at high energies in those 
regions where we cannot extrapolate it. 

The behavior of the extrapolated gaps is essentially lin- 
ear approaching the two second-order phase transitions 
from the outside. This is consistent with a critical expo- 
nent v = 1, characteristic of the two-dimensional Ising 
universality clas o 35 i 36 predicted above. 



B. Phase diagram 

Now we apply the results of the preceding subsection 
to determine the full phase diagram for < J x < J, 
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FIG. 5: (Color online) Numerical results for the phase dia- 
gram in the J2 = J/2 plane. Rung singlet, Haldane, colum- 
nar dimer and staggered dimer phases can be distinguished. 
Crosses have been determined by the phenomenological renor- 
malization group condition 12511 with TVi = 14, TV2 = 16, the 
corresponding lines with TVi = 12, TV2 = 14. For further 
details and the nature of the phase transitions compare the 
text. 



< J± in the plane J 2 = J/2. Numerical scans were 
performed for TV < 14 rungs since they require diagonal- 
izations for a large set of data points (we have performed 
computations for approximately 4000 different pairs of 
J x and J± at TV = 14). In selected regions we performed 
additional computations for TV = 16. Fig. [5] summarizes 
our results. One can distinguish four different gapped 
phases. A rung-singlet and a Haldane phase appear at 
large J± and J x , respectively. Previous results^ for the 
point J x = 0, J± = J, J 2 = J/2 are consistent with 
our conclusion that this lies well inside the rung-singlet 
phase. The rung-singlet and the Haldane phases have a 
non-degenerate ground state for a system with periodic 
boundary conditions. Between these two phases there ex- 
ist two dimerized phases for small values of J x and J±, 
namely a staggered and a columnar dimer phase. The lat- 
ter phases have a two-fold degenerate ground state with 
a spontaneously broken translational symmetry in the 
thermodynamic limit TV — > 00. 

The transition between the staggered and columnar 
dimer phase is of first order and was determined as de- 
scribed in subsection IIV Al level crossings were first de- 
termined for fixed 6 < TV < 14 (16 close to PI) and then 
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FIG. 6: (Color online) Scaled energies of the k = n, k y = 
excited level for Jj_ = 0.9 J, Ji = J/2 as a function of J x . 
Lines connect actual numerical results which are shown by 
symbols. According to the phenomenological renormalization 
group condition 1251 1. crossings of levels for different N are 
estimates for the boundaries of the columnar dimer phase. 

extrapolated to N — > oo using the VBS algorithm. The 
error of this extrapolation is at most on the order of the 
width of the line. This transition line passes through the 
origin, as is expected by consideration of two decoupled 
dimerized chains. 

The transitions from the dimerized phases to the rung- 
singlet or Haldane phases were determined by the phe- 
nomenological renormalization group method l|25() with 
N\ = 12, N 2 = 14 (lines) and N r = 14, N 2 = 16 
(crosses). This is certainly justified for the transition 
between the columnar dimer and rung-singlet phase as 
well as the transition between the staggered dimer and 
Haldane phase, since these should be second-order transi- 
tions belonging to the two-dimensional Ising universality 
class. The transition between the columnar dimer and 
Haldane phase was also estimated by the phenomenolog- 
ical renormalization group method, although this may 
actually be a first-order transition. Nevertheless, finite- 
size effects are small (see e.g. right side of Fig- - There- 
fore, our estimates for the transition between the colum- 
nar dimer and Haldane phase are probably quite accu- 
rate, even if we may not have used the most appropriate 
method. 

The intermediate dimerized phases disappear (at least 
according to the N\ — 12, N 2 — 14 estimates) for large 
interchain couplings at points PI and P2, respectively. 
Beyond the point P2, one expects a direct first-order 
transition from the Haldane to the rung-singlet phase, 
like in the case J 2 = 0i±2As inRef.Q, we determine the 
location of this first-order transition from a cusp in the 
ground-state energy, using the data for N — 14 rungs. 
The transitions between the rung-singlet and Haldane 
phase at J 2 = and J 2 = J/2 probably belong to the 
same universal surface of first-order transitions. 

Finite-size effects of the transition lines are small in 



most regions of Fig. [31 as can be seen by comparison 
of the results for N x = 12, N 2 = 14 and N x = 14, 
N 2 = 16. This applies in particular to the neighbor- 
hood of PI whose location can therefore be considered 
accurate. However, close to P2, we observe finite-size ef- 
fects of the boundaries of the columnar dimer phase. In 
order to discuss this in more detail, Fig. H3 shows scaled 
energies of the k = ir, k y = excited level for J± = 0.9 J, 
J 2 = J/2 and 6 < N < 16. The crossings between the 
scaled levels for different A" estimate the phase bound- 
aries according to the phenomenological renormalization 
group condition (|25|l . Unfortunately there is no system- 
atic dependence of the crossings points on Ni, N 2 such 
that an extrapolation is not possible. Most crossings 
for the transition to the rung singlet phase (left side of 
Fig. |SJ) lie between the Ni = 14, N 2 = 16 crossing at 
J x w 0.642 J and the N ± = 12, N 2 = 14 crossing at 
J x ~ 0.704 J. Furthermore, linear extrapolation of the 
levels for fixed A" suggests a closing of the gap and thus 
a phase transition around J x ~ 0.66 J. Overall, it seems 
likely that the N x — 12, N 2 = 14 data underestimates 
the stability region of the columnar dimer phase, while 
Ni = 14, N 2 = 16 may overestimate it. Thus, the lo- 
cation of the termination point P2 is probably at larger 
J x , J± in the thermodynamic limit than the estimate 
in Fig. Still one can argue that the columnar dimer 
phase has to terminate. Firstly, it cannot extend beyond 
J x = J, since the duality transformation mentioned at 
the beginning of section II I II maps the columnar dimer 
phase which we observe for J x < J to a different phase 
with dimers on the diagonals for J' x = J > J' = J x , 
J' ± = J±. Secondly, the region where J± is much bigger 
than all other exchange constants obviously belongs to 
the rung-singlet phased Therefore, the columnar dimer 
phase has to terminate cither at a point P2 located in the 
region J x < J, or on the line J x = J. If the latter sit- 
uation is realized we cannot exclude the possibility that 
the rung-singlet and Haldane phase are separated by a 
Gaussian quantum critical point located at J x = J. 

Finally, we discuss the relation of the numerical phase 
diagram to the picture obtained by the mapping to the 
Gross-Neveu model presented in subsection IIII Bl Nu- 
merically we see that the columnar dimer phase extends 
to stronger interchain couplings than the staggered dimer 
phase. This gives further support to the observation that 
additional interactions such as ring exchange are nec- 
essary to fine-tune all relevant interchain couplings to 
zero. Indeed, it was found that ring exchange stabilizes 
staggered dimer order ^SSi Therefore, ring exchange will 
cause the points PI and P2 in Fig. [S] to move towards 
each other. PI and P2 will merge when all relevant in- 
terchain couplings are fine-tuned to zero. For the latter 
situation the effective theory of the multi-critical point 
will be described by one gapless Majorana fermion de- 
coupled from the rest of the gapped 0(3) Gross-Neveu 
model as presented in subsection IIII Bl Since the effec- 
tive field theory describing long wavelength fluctuations 
of our model has a form similar to the one considered 



9 



Uj 1.5 

< 



I \ 




1 1 1 


\ \ 

_\ 

\ \ 


— s=o 

-- S=1 
■-■ S=2 




■ \ N 
\ 

~ \ 






\ 

\ ~ ' 
\ 

\ 

\ 

N / 




' ' I ' ' 


i iii 


, , 1 



I I I I I I I l_l I I I I I I I I I 

'.1 1 10 

J, IJ 



FIG. 7: (Color online) The lowest excited levels with total 
spin S = 0, 1, and 2 on a ladder with N = 12 rungs for 
Jx = 0, J2 = J/2 as a function of J± on a semi-logarithmic 
scale. 



been observed for ladders with J2 = J x = (see e.g. Ref. 
Elf). However, in this case none of the bound states has 
lower energy than the lowest edge of the two-triplet con- 
tinuum. Two things happen when one turns on J2 or J x . 
Firstly, the minimum of the single-triplet dispersion gen- 
erally shifts to incommensurate wave vectors. Secondly, 
now a bound state can appear at strong coupling which 
has an energy below the lowest edge of the two-particle 
continuum (see also Ref. ^] for J2 = 0). However, it is 
necessary to go to high orders 4 ^ to obtain quantitative 
agreement with numerical data. Therefore we do not 
pursue this further here. 

Finally, we would like to mention that the bound-state 
nature of the lowest excitation leads to unusual finite-size 
effects close to the phase transitions and thus explains 
the difficulties with the extrapolations encountered in the 
previous subsections. 



V. CONCLUSIONS 



in Ref. 0, the same conclusion concerning the need to 
include additional four-spin interactions in order to sup- 
press all relevant interchain couplings applies to the case 
J2 = considered in Ref. IfJ 



C. Softening of a bound state 

In this last subsection we will comment on the na- 
ture of the elementary excitations. For this purpose we 
start from strong Jj_ where one naturally obtains a rung- 
singlet phased and approach the transition to the colum- 
nar dimer phase along the line J x = 0. Since we just wish 
to make a qualitative point, we restrict to N — 12 rungs 
in this section. Fig.[7|shows numerical results at J x = 0, 
J2 = J/2 for the lowest excited levels in the sectors with 
total spin S — 0, 1, and 2. At strong J± the lowest exci- 
tation is a propagating triplet around AE = Jj_ above a 
non-degenerate ground state. On the other hand, the 
singlet excitation which softens at the transition to a 
columnar dimer state at J± » 0.29 J can be traced to 
the lower boundary the S — 2 continuum which starts at 
AE = 2 J± for J± — > 00. This implies that the singlet 
which at small J± is a low-energy excitation originates 
from a bound state of two triplets at strong J± . For the 
parameters of Fig. [7| this bound state has lower energy 
than the fundamental triplet for J± < 0.8 J and finally 
softens at the Ising transition while states with non-zero 
total spin remain at a finite energy. Similarly, the tran- 
sition from the Haldane phase to the staggered dimer 
phase should occur via softening of an S = two-magnon 
bound state while a single magnon remains gapped at the 
transition. 

It is straightforward to analyze the one-triplet and the 
bound-state problem of two triplets by perturbation the- 
ory at strong J±. In fact, formation of bound states has 



We have studied the ground-state phase diagram and 
the nature of elementary excitations of a frustrated spin 
ladder. The effect of relevant and marginal interchain 
exchanges has been analyzed by an effective field theory 
approach. Particular attention was paid to the interplay 
between two clearly separate scales set by inchain frustra- 
tion and interchain couplings. First, a rung-singlet and a 
Haldane phase are known to appear for J2 = Q^Li^i^JL 
We predicted that for sufficiently large J2, two dimcr- 
ized phases appear between the rung-singlet and Hal- 
dane phase. The transitions from the dimerized phases 
to the rung-singlet and Haldane phases are found to be 
in the two-dimensional Ising universality class, whereas 
the transition between the columnar and staggered dimer 
phases is of first order. 

Furthermore, we have verified these predictions by ex- 
act diagonalization. The full phase diagram has been 
determined in the plane J2 = J/2, see Fig. [5J The 
staggered dimer phase was found to terminate first at 
a point PI while the columnar dimer phase extends to 
stronger interchain coupling, terminating at P2. Beyond 
P2 there is a direct first-order transition from the rung- 
singlet to the Haldane phase if P2 is located at J x < J, 
like for J 2 = 0i±2 Should P2 be located at J x = J, 43 
rung-singlet and Haldane phase might be separated by a 
second-order phase transition point. 

The part of the phase diagram with J x > J is obtained 
from the one with J x < J by a duality- type argument. 13 
In particular, the columnar and staggered dimer phases 
map to two different dimerized phases with dimers on the 
diagonals. 

Finally, we have observed an interesting phenomenon 
when approaching the columnar or staggered dimer phase 
from the rung-singlet or Haldane phase, respectively. 
Namely, a singlet bound state forms below the two-triplet 
scattering continuum. This singlet bound state crosses 
the fundamental triplet excitation and finally softens at 
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the Ising phase transition. 

Since all phases are gapped, they will be stable un- 
der weak higher-dimensional coupling between the lad- 
ders. If one approaches the J\ — J2 — J3 frustrated square 
lattic o 9 i 10 by such inter-ladder coupling, the rung singlet 
phase will be adjacent to a Neel phase, and the Haldane 
phase will share a boundary with a phase where spins 
are ordered antiferromagnetically along one 'ladder' di- 
rection and ferromagnetically along the other direction. 
Our results show that in such a generalized phase di- 
agram there will be further dimerized phases, in par- 
ticular some with two-dimensional columnar and stag- 
gered dimer patterns. We therefore hope that our results 
will also be useful to elucidate the phase diagram of the 
J 1 — J-2 — J3 frustrated spin- 1/2 Heisenberg antiferromag- 
net on the square lattice 
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